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Abstract 

In this paper, we prove some general convergence theorems for the Picard 
iteration in cone metric spaces over a solid vector space. As an application, we 
provide a detailed convergence analysis of the Weierstrass iterative method 
for computing all zeros of a polynomial simultaneously. These results improve 
and generalize existing ones in the literature. 
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1. Introduction 

In the first part of the paper, we study the convergence of the iterative 
processes of the type 


Xn+l=TXn, n = 0,1,2,..., 


( 1 . 1 ) 


where T: H C A ^ A is an iteration function in a cone metric space (A, d) 
over a solid vector space (V, Y). Cone metric spaces have a lon g hi story (see 
Collatz j^, Zabrejko ji^, Jankovic, Kadelburg and Radenovic 10|, Proinov 
29j and references therein). For an overview of the theory of cone metric 
spaces over a solid vector space, we refer the reader to and j^. Section 2]. 

In the second part of the paper, we study the convergence of the famous 


Weierstrass method j39| for computing all zeros of a polynomial simultane¬ 


ously. This method was introduced and studied for the hrst time by Weier¬ 
strass in 1891. In 1960-1966, the method was rediscovered by Durand (in 
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implicit form), Dochev Kerner [llj] and Presic j23|]. For this reason, it is 


also known as ‘Durand-Kerner method’, ‘Weierstrass-Dochev method’, etc. 
For an overview of iterative methods for simultaneous hnding of polynomial 
zeros, we refer the reader to (lil . 14, 18, 34 • 


Throughout this paper, (K, I'D denotes an arbitrary normed (valued) 
held with absolute value | • |, and ]K[z] denotes denotes the ring of polynomials 
in one variable 2 ; over IK. Let / G IK)^;] be a polynomial of degree n >2. We 
consider the zeros of / as a vector in K”. More precisely, a vector ^ G K" is 
said to be a root-vector of / if 


/(^) ^ aoP[(^ - ^i) for all ^ e K, 


( 1 . 2 ) 


2=1 


where Oq G K. Obviously, / has a root-vector in K" if and only if / splits in 
K. Recall that the Weierstrass method is dehned by the following iteration 

-W{x^), A: = 0,1,2,..., (1.3) 


where hF: "D C K” —)■ K” is dehned by W{x) = (IFi(x),..., Wn{x)) with 

-- (r = l,...,n), (1.4) 

ao J_J_(a(* - Xj) 
j¥=i 

where Oq is the leading coefficient of / and "D is the set of all vectors in 
K” with distinct components. The operator W is called the Weierstrass 
correction. Sometimes we write Wf instead of W to indicate that the operator 
W is generated by /. It is easy to see that the Weierstrass correction Wj is 
invariant with respect to multiplication of / by a non-zero constant c G K. 
Obviously, the Weierstrass iteration fll.3p can be represented in the form fll.ll) 
with the iteration function T: T> G K"' —)■ IK"" dehned by 


T{x) = X — W{x). 


(1.5) 


The aim of this paper is twofold. First, we present some general con¬ 
vergence theorems with error estimates for the Picard iteration fll.ll) . These 


results extend some of the results in [27|, 12^. Second, using these results we 


provide a detailed convergence analysis of the Weierstrass method fll.ip . The 
new results for the Weierstrass method improve the corresponding results of 


0, i, 0 0 0, EsH 0, S 10 s 0, Bie 
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The paper is structured as follows: 

In Section [21 we present some preliminaries results and notations that 
will be useful in the sequel. 

In Section [3l we establish two general convergence theorems with error 
estimates for iterated contractions at a point in cone metric spaces. The hrst 
one extends Theorem 3.6 of 27 


In Section 01 we establish two general semilocal convergence theorems 
with error estimates for iterative processes of the type fll.ll) . These results 
extend Theorems 5.4 and 5.6 of [2^. As a consequence we obtain a conver¬ 
gence theorem with error estimates for iterated contractions in cone metric 
spaces, which extends Theorem 6.5 of j^. All results in this sections are gen¬ 
eralizations of the Banach Contraction Principle [ij as well as of the Iterated 
Contraction Principle given in [l^, Chap. 12] and j^. 

In Section [5l we present some inequalities in K"' and notations which will 
be useful in the next sections. 

In Section [6l we obtain a local convergence theorem with error esti¬ 
mates for the Weierstrass method which improves the results of Dochev |4|, 
Kyurkchiev and Markov 13|, Yakoubsohn ^ and Proinov and Petkova [32 1. 

In Section [71 we obtain another local convergence theorem with error 
estimates for the Weierstrass method which improves and generalizes the 
results of Wang and Zhao 37], Tilli 36| and Han [^. 

In Section [HI we prove a new convergence theorem for the Weierstrass 
method under computationally verihable initial conditions. The main result 
of this section generalizes, improves and complements all previous results 
in this area, which are due to Presic j^, Zheng (iR 4^ . Wang and Zhao 
38l . [441 . Petkovic, Carstensen and Trajkovic ll9|, Petkovic 


__ Petkovic, 

Herceg and Ilic 211, Batra j^, Han j^, Petkovic and Herceg |20| and Proinov 


24j |. The new result in this section also gives computationally verihable error 


estimates, a localization formula for the Weierstrass iteration (in as well 
as a sufficient condition for a polynomial / G K[ 2 ] of degree n > 2 to have n 
simple zeros in the held K. Finally, we provide an example which shows the 
exactness of the error estimates of our semilocal theorem for the Weierstrass 
iterative method. 


2. Preliminaries 

Throughout the paper, J denotes an interval in M+ containing 0, that is, 
an interval of the form [0,i?], [0, i?) or [0,cxo), where i? > 0 . We use the 
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abbreviation for the nth iterate of a function ip\ J ^ J. For n G N we 
denote by Sn(t) the following polynomial 


n—1 

S„(i) = E*‘- (2-1) 

k=0 

If the case n = 0, we set S'o(t) = 0. Throughout the paper, we assume by 
dehnition that 0 ° = 1 . 


Definition 2.1 ( 2^). A function 9 ?: J —> M+ is called quasi-homogeneous 
of degree r > 0 on J if it satishes the following condition 


< y for all A G [0,1] and t G J. (2.2) 


It is easy to prove that a function : J —)■ M+ is quasi-homogeneous 
of degree r > 0 on J if and only if there exists a nondecreasing function 
$ : J — )■ M+ such that (p(t) = ^{t) for all t G J. 

Let us give an example for quasi-homogeneous functions, which we will 
use in Sections E] and [71 

Example 2.2. Let n G N and be a quasi-homogeneous function of degree 
r > 0 on an interval J. Then the function $ dehned by 


is also quasi-homogeneous of degree r on J. 

Proof. It follows from the identity = J2k=i because the sum 

of quasi-homogeneous functions of degree r on J is quasi-homogeneous of 
degree r on J. □ 

Definition 2.3 ([^). A function 9 ?: J —)■ J is called a gauge function of 
order r > 1 on an interval J if it is quasi-homogeneous of degree r on J and 


ip{t) < t for all t E J. (2.3) 

A gauge function ip of order r on J is said to be a strict gauge function if 
the inequality in fl2.3p holds strictly whenever t > 0. 

Proposition 2.4 ([^). If<p'. J K+ is quasi-homogeneous of degree r > 1 
on an interval J and R > 0 is fixed point ofip in J, then ip is a gauge function 
of order r on [0, R], Moreover, if r > 1, then ip is a strict gauge function of 
order r on [0, R). 
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Definition 2.5 (|^). Let T: D C X ^ X he a map of a set X. A function 
E: D ^ M+ is said to be function of initial conditions of T (with a gauge 
function ip on J) if there exists a function p: J ^ J such that 

E{Tx) < p{E{x)) for all x G -D with Tx ^ D and E{x) G J. 

Definition 2.6 ([^). Let T: D G X ^Xbea map of a set X and 
E : D ^ M+ be a function of initial conditions of T with a gauge function 
on J. Then a point x E D is said to be an initial point of T ii E{x) E J and 
T”x G D for all n > 0. 

Proposition 2.7 (j^). Let T: D G X ^ X be a map of a set X and 
E: D M+ be a funetion of initial conditions of T with a gauge function p 
on J. Suppose 


X E D with E{x) G J implies Tx E D. (2.4) 


Then every point xq E D such that E{xq) E J is an initial point ofT. 

Proposition 2.8 f(^L Let T D G X ^ X be a map of a set X and 
E: D ^ M_i_ be a function of initial conditions of T with a gauge function 
p on an interval J. If x E D is an initial point of T , then every iterate 
Xn = T”x {n = 0,1,2,...) is an initial point of T. Besides, if p is a gauge 
function of order r >1, then for all n > 0, 

Eixn+i) < X"E{xn) and E{xn) < E{xo) 

where A = 0(£'(xo)) and 0: J —)■ [0,1] zs a nondecreasing function such that 


pit) = f0(f) for all t E J. 


(2.5) 


The following definition extends Definition 3.4 of j28|. For an ordered 


vector space {Y, X) we denote by Y+ the positive cone of Y, that is, 


Y+ = {x E Y : X ^ 0}. 


Definition 2.9. Let {X,d) be a cone metric space over an ordered vector 
space (y, y). Suppose T: D G X ^ X is an operator and E: D ^ M+ is 
a function of initial conditions of T with gauge function p on an interval 
J. Then a function E: D ^ Y+ is said to be a convergence function of T 
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(with respect to E and with control functions (3 and 7) if there exist two 
nondecreasing functions /9: J —?■ [0,1) and 7: J —)■ M+ such that 

F{Tx) :< (3{E{x)) E{x) for all x E D with Tx E D and E{x) E J, (2.6) 

d{x,Tx) :< 'j{E{x)) E{x) for all x E D with E{x) E J. (2.7) 

Usually we need stronger conditions for the control function /3. Namely, 
we assume that /3: J —)■ [0,1) is a nondecreasing function satisfying the fol¬ 
lowing two conditions: 

t(3{t) is a strict gauge function of order r on J, (2.8) 

WtEJ: = 0 implies (3{t) = 0, (2.9) 

where 0: J ^ [0,1] is a nondecreasing function satisfying fl2.5p . 

It is easy to see that condition fl2.9p is equivalent to the existence of a 
function 0: J —)■ M+ such that 


(3{t) = for all t E J. 


The following proposition extends Lemmas 3.5 and 3.8 of 28 


( 2 . 10 ) 


Proposition 2.10. Let T: D G X ^ X be an operator of a cone metric 
space {X, d) over an ordered vector space {Y,^),E:D^ M_|_ be a function 
of initial conditions ofT with gauge function ip on J satisfying fl2.3p . and let 
F; H —)■ y_|_ be a convergence function of T with control functions (3 and 7. 
If X E D is an initial point of T, then for all n > 0, 


E{xn+i) Y hE{xn) and E{xn) Y h'^E{xo), ( 2 - 11 ) 

where Xn = T'^x and h = (3{E{xo)). Besides, if (p is a gauge function of 
order r > 1 on J and (3 satisfies fl2.8p and fl2.9l) . then for all n >0, 

E{xn+i) Y eX^"E{xn) and E{xn) Y e^X^-(^^E{xo), (2.12) 

where X = (j){E{xo)), 6 = 'f{E{xo)) and if: J ^ R+ is a function satisfying 

dnnp. 
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Proof. We shall prove only fl2.12p since the proof of fl2.1ip is similar. Suppose 
Xo is an initial point of T. We have Xn & D and E{xn) € J since every iterate 
Xn is an initial point of T. Setting x = x„ in fl2.6p . we obtain 

F{xn+i) ^ l3{E{xn)) E{xn). (2.13) 

Condition fl2.8l) implies that /3 is a quasi-homogeneous function of degree 
r — 1 on the interval J. From this and Proposition 12.81 we get 

P(E(x„)) < = ey". 

(2.14) 

Combining fl2.13p and fl2.14p , we get the hrst inequality in fl2.12p . The second 
inequality in fl2.12p follows from the hrst one by induction on n > 0. □ 

We end this section with two proposition from the theory of cone metric 
spaces over a solid vector space. These statements are trivial in the case of 
metric spaces. 

Proposition 2.11 f[^). Let {X, d) he a cone metric space over a solid 
vector space (Y, Y)- Suppose (xn) is a sequence in X satisfying 

d{xn,x) P bn for all n > 0, (2-15) 

where x is a point in X and {bn) is a sequence in Y converging to 0. Then 
the sequence (x„) converges to x. 

Proposition 2.12 (j^). Let {X,d) he a cone metric space over a solid 
vector space {Y, Y)- Suppose (x^) is a sequence in X satisfying 

d{xn,Xm) Y bn for all n,m>0 with m > n, (2-16) 

where {bn) is a sequence in Y which converges to 0. Then: 

(i) The sequence (x„) is a Cauchy sequence in X. 

(ii) If {xn) converges to a point x E X, then 

d{xn,x) Y bn for all n>0. (2-17) 
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3. Local convergence theorems in cone metric spaces 


In this section, we present two general convergence theorems with error 
estimates for iterated contractions at a point in cone metric spaces. 

Definition 3.1 (Iterated Contraction at a Point). Let T: D G X ^ X he 
an operator of a cone metric space {X,d) over a solid vector space (Y, Y), 
and let E: D ^ M_|_ be a function of initial conditions of T with a gauge 
function on an interval J. Then T is said to be an iterated contraetion with 
respect to E at a point ^ E D (with control function (3) if E{^) G J and 

d{Tx,^) Y P{E{x)) d{x,^) for all x G D with E{x) G J, (3.1) 

where /3: J ^ [0,1) is a nondecreasing function. 

Proposition 3.2. LetT: D G X ^ X be an operator of a cone metric space 
{X, d) over a solid veetor space (L, ^), and let E: D —)■ M_|_ be a function of 
initial conditions of T with a gauge function on an interval J. If T is an 
iterated contraction with respect to E at a point f & D, then ^ is a unigue 
fixed point ofT in the set U = {x E D : E{x) G J}. 

Proof. Setting x = .^ in 03.11) . we get d{T^,^) P 0. Therefore, d{T^,^) = 0 
which means that is a fixed point of T. Suppose t] E U is also a fixed point 
of T. Applying 03.11) with x = rj we obtain d{ri,^) Y dijj^ff) which 

implies that d{ri,^) Y 0 because the values of (3 are less than 1. From this, we 
conclude that rj = ^. Hence, ^ is a unique fixed point of T in the set U. □ 


In the following two theorems, we consider the problem of approximating 
the fixed points of iterated contractions at a point in a cone metric space. 
The first one extends Theorem 3.6 of 


27 


Theorem 3.3. Let T: D G X ^ X be an operator of a cone metric space 
{X,d) over a solid vector space (T, ^), and let E: D —)■ M_|_ be a function 
of initial conditions of T with a gauge function g) of order r on an interval 
J. Suppose T is an iterated contraction with respect to E at a point with 
control function f3 satisfying 02.Sp and 02.9p . Then for each initial point xq 
ofT the Picard iteration 01.ip remains in in the setU = {x E D : E{x) E J} 
and converges to f with error estimates 

d(x„+i,0 ^ 6'A'’"d(x„,0 and d(x„,0 ^ 0 (3-2) 


for all n > 0, where A = (j){E{xo)), 9 = fj^E^xo)) and J ^ M+ is a func¬ 
tion satisfying 02 .10p . 















Proof. It follows from the second estimate in fl3.2p and Proposition 12. Ill with 

that {xn) converges to ^ in X. By the dehnition of initial points, it follows 
that Xo G U. Starting from xq the iterative sequence fll.ll) remains in the set 
U because every iterate is an initial point of T. Setting x = x„ in fl3.ip . 
we obtain 

d{Xn+l,0 ^ f^iEiXn)) d{Xn, 0 , (3-3) 

From 03.31) and 02.141) . we get the hrst estimate in 03.21) . The second estimate 
in 03.21) is a consequence of the hrst one. □ 

Setting = 0 in Theorem 13.31 we get the following result. 

Corollary 3.4. Let T: D G X ^ X be an operator of a cone metric space 

{X, d) over a solid vector space (F, ^), and let E: D ^ M+ be a function of 

initial conditions ofT with a strict gauge function (p of order r on an interval 
J. IfT is an iterated contraction with respect to E at a point ^ with control 
function (f satisfying 02.51) . then for each initial point xq of T the Picard 
iteration 01.11) remains in in the set U = {x E D : E{x) E J} and converges 
to ^ with error estimates 

d{Xn+l,0 ^ d{Xn,0 d{Xn,0 ^ d{xo,0 (3-4) 

for all n>0, where A = (j){E{xQ)). 

Theorem 3.5. Let T: D G X ^ X be an operator of a cone metric space 

{X, d) over a solid vector space {Y, and let E: D M_|_ be a function of 

initial conditions of T with a nondecreasing gauge function p on an interval 
J. Suppose T is an iterated contraction with respect to E at a point ^ with 
control function fd. Assume there exist a E (0,1), r > 1 and a nondecreasing 
function c: [0, a] —?■ J such that 

/9(c(t)) < t and p{c{t)) < c{f) for all t E [0,cr]. (3.5) 

If xo is an initial point point of T satisfying 

E{xo) < c{a), (3.6) 

then Picard sequence 01.ip converges to ^ with error estimates 

d{xn+i,i) < d{xn,i) and d(xn, ^ ^ c?(a^o, 0 (3-7) 
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Proof. We shall prove only the hrst estimate in fl3.7p . Using Dehnition 12.51 
and the second inequality in fld.Sp it is easy to see that for every initial point 
X oi T and t e [0, a], 

E{x) < c{t) implies E{Tx) < c{f). 

Then by induction one can prove that 


E{xn) < c (cr'’ ) for all n >0. 

This inequality together with the hrst inequality in fl3.5p implies 

/3{E{Xn)) < 

From this and fl3.3p . we get the hrst estimate in fl3.7p . □ 


4. General convergence theorems in cone metric spaces 


In this section, we establish two general convergence theorems with er¬ 
ror estimates for iterative processes of the type (II.ip . These results extend 


Theorems 5.4 and 5.6 of 28 


Theorem 4.1. LetT: D G X ^ X be an operator of a complete cone metric 
space {X,d) over a solid vector space (U, E: D ^ M_|_ be a function of 
initial conditions ofT with gauge function if on an interval J satisfying fl2.3p . 
and let E: D ^ Yj^ be a convergence function of T with control functions f3 
and 7. Then the following statements hold true: 

(i) Convergence. Starting from any initial point xq of T, the Picard 
iteration fll.ip is well-defined, remains in the closed ball U{xo,p) and 
converges to a point ^ G U{xq, p), where 


P 


l{E{xo)) 

l-(5{E{xo)) 


F(xo). 


(ii) A PRIORI ESTIMATE. For all n >0 we have the following estimate 


d{xn, 0 ^ 


where h = {3{E{xo)). 
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(iii) First a posteriori estimate. For all n > 0 we have the following 
estimate 

(iv) Second a posteriori estimate. For alln > 0 we have the following 
estimate 


(v) Existence of a fixed point. If f, e D and T is continuous at 
then f, is a fixed point ofT. 

Proof. It follows from fl2.3p and Proposition 12.81 that E{xn) < E{xq) for all 
n > 0. Let m, n G M with m > n. From the triangle inequality, fl2.7p and 
Proposition 12.101 we get 

mm m 

d{Xn,Xm) '^d{Xj,Xj+i) F 

j=n j=n j=n 

x E'‘d l(SK))F(io) S y3^7(B(x„))F(x„). 

\j=n / 

Therefore, 


d{xn,Xm)^bn, where bn = - -r 7(-^(^o))-^(a^o)- (4.4) 

1 — h 

Note that (6„) is a sequence in Y which converges to 0 since h” ^ 0 in M. 
By Proposition 12.12f ii. we conclude that (x^) is a Cauchy sequence. By the 
completeness of the space X, we deduce that (x^) converges to a point ^ E X. 
Now it follows from Proposition I2.12f iii that 


d'i.Xny 0 — bn 


(4.5) 


for every n > 0. From this we get the estimate fl4.ip . Setting n 
we get 


d{xo,^) Y 


7{E{xo)) 

l-(3{E{xo)) 


E{xo) 


0 in dH]), 
(4.6) 
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which means that ^ G U{xo,p). The inequality fl4.6p holds for every initial 
point Xq of T. Therefore, applying fl4.6p to x^ instead of Xq, we obtain fl4.2p . 
From fl4.2p and fl2.6p . we obtain 


d{xn+i,i) ^ 


^{E{Xn+l)) 


l-P{E{Xn+l)) 

(3{E{Xn)) 


l-P{E{xn+i)) 


E{Xn+l) 

'y{E{Xn+l)) E{Xn) 


which proves the estimate (14.31) . Setting n = 0 in (14.4p . we get d{xo, Xm) ^ P 
for every m > 0. Hence, the sequence {xn) lies in the ball U{xo,p). □ 


Theorem 4.2. Let {X, d) be a complete cone metric space over a solid vector 
space {Y, , T: D G X ^ X be an operator, E: D ^ M+ be a function of 

initial conditions ofT with a gauge function p of order r >1 on an interval 
J, and let E ■. D Y+ be a convergence function of T with control functions 
(3 and'y satisfying (12.Sp and (12.9p . Then the following statements hold true: 

(i) Convergence. Starting from any initial point xq of T, the Picard 
iteration (II.ip is well-defined, remains in the closed ball U{xo,p) and 
converges to a point f G U{xo,p), where 


P 


imxo)) 

l-l3{E{xo)) 


F{xo). 


(ii) A PRIORI ESTIMATE. For alln > 0 we have the following error estimate 



(iv) 


an\Sn(r) 

d(x„,C S FM. U7) 

where A = f^E^xo)), 9 = fj^E^xo)), 0: J —)■ [0,1] is a nondecreasing 
function satisfying (12.Sp and fj-. J ^ M+ is a function satisfying (12.101) . 
First a posteriori estimate. For all n > 0 we have the following 
error estimate 


d{xn,0 ^ 


7{E{xn)) 

l-(3{E{Xn)) 


E{Xn). 


(4.8) 


Second a posteriori estimate. For alln > 0 we have the following 
error estimate 


d{Xn+l, 0 ^ 1 ]r 7(E(Xn+l)) F(Xn). (4.9) 


where = (j){E{xn)) and 9^ = 'f{E{xn)). 
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(v) Some other estimates. For all n > 0 we have 

F{xn+i) < e F{x^) and F(x„) ^ F(xo). (4.10) 

(vi) Existence of a fixed point. If ^ e D and T is continuous at f, 
then f is a fixed point ofT. 

Proof. Conclusions (i), (iii) and (vi) follow immediately from Theorem 14.11 
Conclusion (ii) follows from fl4.8p . Proposition 12.81 Proposition 12. lUl and in¬ 
equality fl2.14p . Conclusion (v) follows from Proposition 12.101 It remains to 
prove (iv). It follows from fl2.10p that 

l3{E{Xn)) = e^Xn. (4.11) 


On the other hand, taking into account that E{xn+i) < <p(E(a;„)) and that 
(3 is quasi-homogeneous of degree r — 1 on J, we get 


(3{E{Xn+l) < l3{p{E{Xn)) = (3{XnE{Xn)) < (A 
Now conclusion (iv) follows from 04.31) . 04.lip and 04.12p . 


nr^l3{E{Xn)) = 0n(An)C 

(4.12) 
□ 


Definition 4.3 (Iterated Contraction 28|]). Let T: D C X —)■ X be an op¬ 


erator of a cone metric space (X, d) over a solid vector space (Y, Y), and let 
E: D ^ M_i_ be a function of initial conditions of T with a gauge function on 
an interval J. Then T is said to be an iterated contraction with respect to E 
(with control function fi) if 


d{Tx,T‘^x) Y j3[E{x)) d{x,Tx) for all x G D with x G D and E{x) G J, 


where /3: J ^ [0,1) is a nondecreasing function. 


Setting E{x) = d{x,Tx) in Theorem l4.2l we get the followin g co nvergence 


result for iterated contractions, which extends Theorem 6.5 of [28 


Corollary 4.4. Let T: D G X ^ X be an operator of a cone metric space 
(X, d) over a solid vector space (X, ^), and let E: D —)■ M_|_ be a function 
of initial conditions of T with a gauge function ip of order r on an interval 
J. If T is an iterated contraction with respect to E with control function fi 
satisfying 02.81) and 02.91) . then the following statements hold true: 
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(i) Convergence. Starting from any initial point xq of T, the Picard 
iteration fll.ip is well-defined, remains in the closed ball U{xo,p) and 
converges to a point f G U{xo,p), where 

_ d{xQ,Txo) 

^ l-/3{E{xo))' 

(ii) A PRIORI ESTIMATE. For alln> 0 we have the following error estimate 

Qn^Sn(r) 

d{xn,i) A d{xo,Txo), 



where A = ^{E^xq)), 6 = fi^E^xo)), 0: J —)■ [0,1] is a nondecreasing 
function satisfying fl2.5p and J ^ M+ is a function satisfying (12.101) . 
First a posteriori estimate. For all n > 0 we have the following 
error estimate 


d{,Xn,i) A 


d(Xn, 3^n+l) 

l-/3(F(a;J)- 


(iv) Second a posteriori estimate. For alln > 0 we have the following 
error estimate 


d{xn+i,i) A 


9n\ 


n 


1 - OniXnY 


d{Xn, Xn-\-Y ■ 


where A„ = (j){E{xn)) and 9n = 'fi{E{xn))- 
(v) Some other estimates. For all n>0 we have 

d{xn+i,Xn+ 2 ) A 9 X'" d{xn, Xn+i) and d{Xn, Xn+l) A 6*” d{xo, Xi). 


(vi) Existence of a fixed point. IfU{xo,p) c D andT is continuous, 
then f is a fixed point ofT. 


Remark 4.5. Each of the theorems and corollaries of this section is a gen¬ 
eralization of the Banach Contraction Principle [i| (see also 291) a s well as 


of the Iterated contraction principle given in [15|, Chap. 12] and [29|. For ex¬ 
ample, Corollary 14.41 with E{x) = d(x, Tx), ip{t) = Xt and (3{t) = A, where 
A G [0,1), yields Theorem 10.1 of 
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5. Some inequalities in 


In this and the next sections, we use the following notations and conven¬ 
tions. The vector space K"’ is equipped with the p-norm 


X 



for some 1 < p < oo. 


(5.1) 


The Banach space (M"', || • ||p) is equipped with the coordinate-wise ordering 
■< dehned by 


X ^ y if and only if Xj < yi for each i E In- (5.2) 

Here and in what follows, we denote by In the set of indices 1,... ,n, i.e. 
In = {I,, n}. It is easy to see that (M", || • ||p, is a solid vector space. 
We define the map || ■ ||: K"' —)■ M" by 


ll^ll = (kl|.'">kn|) 


(5.3) 


Then (K”, || • ||) is a cone normed space over M"'. 

X 

Furthermore, for two vectors x G IK” and p G M” we denote by — a vector 

y 

in M” dehned by 

'|Xi| \Xn 


X 

y 


(5.4) 


Pi Vn 

provided that y has only nonzero components. 

We use the function d\ K” —)■ M” dehned by d{x) = (di(x), ... ,dn{x)) 
with 

di{x) = mm\xi — Xj\ (f = l,...,n). (5.5) 

Proposition 5.1 f(3ol|b Let u,v E K” and 1 < p < oo. If the vector v has 
distinet components, then for all i,j E In the following two inequalities hold: 


(i) \ui-Uj\ > 1 -2^/'' 


u — V 


d{v) 


Vi - V 


3 \’ 


(ii) \ui-Vj\> 1 - 


u — V 


d{v) 


\Vi - Vi 


The next two lemmas are immediate consequences of the previous one. 
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Proposition 5.2. Letu,v G K” and 1 <p < oo. If the vector v has distinct 
components, then 


d{u) ^ 1 - 2 ^/-? 


u — V 


d{v) 


d{v). 


Proposition 5.3. Let G K” and 1 <p < oo. If the vector v has distinct 
components and 


u — V 


d{v) 


< 


1 


21/q ’ 


then the vector u also has distinct components. 


Given —oo < r < oo, we define the power mean function : 


by 


l/r 


Xi 


Mr{x)= - 5 ^ 

\ i=l 

The value of Mr{x) for r = 0, ±oo is assumed to be the limit as r —)■ 0, ±oo. 
Proposition 5.4 (Power Mean Inequality). If —oo <r < s < oo, then 
Mr{x) < Ms{x) for every a: G IK"', 

and the equality holds only if all the components of x are equal to each other. 

The next proposition can easily be proved by the power mean inequalities 
Mq < Ml and Mi < Mp for p > 1. 


Proposition 5.5 ([30, l32j). Let u G K" and 1 < p < oo. Then 


]^(1 + Ui) 


2=1 


< 1 + 


\u\ 


n 


i/p 


and 


]^(1 + Ui) — 1 


2 = 1 


< 1 


\u\ 


n 


i/p 


- 1 . 


Throughout the next sections, for a given 1 < p < oo, we always denote 
by q the conjugate exponent of p, that is. 


1 < g < oo and —I— = 1. 

p q 


(5.6) 
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6. Local convergence of the first kind of the Weierstrass method 


Let / G K[z] be a polynomial of degree n > 2 which has n simple zeros in 
IK and be a root-vector of /. In this section we stndy the convergence of the 
Weierstrass method fll.Sp with respect to the fnnction of initial conditions 
L: K"' —)■ M_|_ dehned by 




x-i 


d(i) 


(6.1) 


where 1 < p < cx). We prove that the Weierstrass iteration fnnction T dehned 
by (II.5p is an iterated contraction with respect to E at the point As a 
resnlt, we obtain a local convergence theorem with error estimates for the 
Weierstrass method, which improves the resnlts of Dochev Kynrkchiev 


and Markov 13 , Yakonbsohn 


and Proinov and Petkova 32 


Lemma 6.1. Let f G K[z] be a polynomial of degree n > 2 whieh splits in 
K, ^ G K” be a root-veetor of f, x & K” and 1 < p < oo. Then for i ^ j, 


Xi — Xj\ > {1 — 2^^'^E{x)) dj{f) and \xi — > {1 — E{x)) di{f), 


where E: K" —)■ M+ is defined by (16. ip . 

Proof. Setting in Proposition 15.11 u = x and v = f and taking into acconnt 
the dehnition of d{f), we obtain the statement of the lemma. □ 

Lemma 6.2. Let f G ]K[z] be a polynomial of degree n >2 which splits in K, 
f be a root-vector of f, and let x G K” be a vector with distinct components. 
Then for every i G In, 


Ti(x) - 



( 6 . 2 ) 


where 


Uj = 


^3 - 

Xi — Xj 


(6.3) 


Proof. Let i G /„ be hxed. Taking into acconnt the identity fll.2l) . we obtain 


Wfx) 


[Xi 



Xi 

Xi Xj 


3¥=i 


From this and (II.5p . we get (16.2p . 


□ 
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Lemma 6.3. Let f G ]K[z] be a polynomial of degree n > 2 which has n 
simple zeros in K, ^ be a root-vector of f and 1 < p < oo. Suppose a vector 
a: G K” satisfies 


E{x) 


x-f 

diO 


1 


(6.4) 


where the function E is defined by fl6.ip . Then x eT>, 

E{Tx) < ip{E{x)) and ||Tx — ^ (/)(L'(a;)) ||a; — (6.5) 


where T is the Weierstrass iteration function defined by fll.5p . and the real 
functions ip and (j) are defined by 

(^+ (^_i)i/p|i_2i/n) ) ip{t)=t(j){t). (6.6) 

Proof. It follows from Proposition 15.31 with u = x and v = f that x E TP. Let 
i G In be hxed. Combining Lemma 16.21 and Proposition 15.51 we obtain 


Ti{x) - < 




(n - ly/p 


n—1 

- 1 


\Xi 


(6.7) 


where u = {uj)j^i G K” ^ and Uj is dehned by (16.31) . It follows from 
Lemma 16.11 that 


\uj\ — 





X 2 X j 

^ {\-2^I^E{x))d^{x) 


( 6 . 8 ) 


which yields 


,11 ^ r{x) 

W - 1 - 21/«£’(i)' 


From (16.7p and the last inequality, we get 


\Ti{x) - fi\ < 4){E{x))\xi - fi\ (6.9) 

which yields the second inequality in (16.5p . The hrst inequality in (16.5p 
follows from (16. 9 p dividing by dff) and taking the p-norm. □ 

Now we are ready to state the main result of this section. 
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Theorem 6.4. Let f G ]K[z] be a polynomial of degree n > 2 which has n 
simple zeros in K, ^ be a root-vector of f and 1 < p < oo. Suppose G K"" 
is an initial guess such that 


E{x^) 


dio 


< R{n,p) 

V 


2i/(”-i) — 1 

21/9 (2i/(—1) - 1) + (n - 1 )-Vp ’ 


( 6 . 10 ) 


where the function E is defined by flO.ip . Then the Weierstrass iteration 
oa is well-defined and converges quadratically to ^ with error estimates 


_ ^11 ^ llx*^ - ell and llx'^-ell ^ A2'“-i||a:°-e|| (6.11) 


for all k >0, where A = fi^E^x^)) and the real function cf) is defined by 06.61) . 

Proof. By Example 12.21 the function if defined by 06.61) is quasi-homogeneous 
of the second degree on [0, l/2i/9). On the other hand, R = R{n,p) is a hxed 
point of f in (0,l/2i/9). Then according to Proposition 12.41 is a strict 
gauge function of the second order on the interval J = [0, R). 

Now we shall apply Corollary 13.41 to the Weierstrass iteration function 
T: V C. —)■ It follows from Lemma 16/^ that: 

• E: P —)■ R_|_ is a functions of initial conditions of T with gauge function 
f of order r = 2 on J. 

• T; P —)■ KL is an iterated contraction with respect to E at ^ with control 
function <p. 

It remains to prove that every vector x^ with E{x^) G J is an initial 
point of T. According to Proposition 12.71 it is sufficient to prove that 


X eV with E{x) G J implies Tx G V. (6.12) 

Let X eR he such that E{x) E J. From x eR, we get Tx E K". From the 
hrst inequality in 06.51) taking into account that E{x) G J and f. J —?■ J, we 
conclude that E{Tx) G J. Then applying Lemma 16/^ to Tx instead of x, we 
deduce Tx eR which proves 06.12p . 

Now the statement of the theorem follows from Corollary 13.41 □ 

Let 0 < h < 1 be a given number. Solving the equation (/)(t) = h in the 
interval (0, R{n,p)) we can formulate Theorem 16.41 in the following equivalent 
form. 
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Theorem 6.5. Let f G ]K[z] be a polynomial of degree n > 2 which has n 
simple zeros m K, ^ a root-vector of f, 1 < p < oo and 0 < h <1. Suppose 
G IK” is an initial guess such that 


diO 


< R{n,p, h) 


(1 + hfSn-l) _ 1 

2V? ((1 + /i)V(n-l) _!) + („_ l)-l/p • 


(6.13) 


Then the Weierstrass iteration fll.3p is well-defined and converges guadrati- 
cally to ^ with error estimates 


\x‘ 


k+l 


— ^\\ ^ h'^ 11^^ “'Cll o,nd \\x^ — ^\\ h'^ 


\x 


-ell (6.14) 


for all k > 0. 

In the next two corollaries we denote by sep(/) the separation number of 
/ which is dehned to be the minimum distance between two zeros of /, that 
is, 

sep(/) = minimiI. (6.15) 

Corollary 6.6. Let f G IK).^] be a polynomial of degree n > 2 which has 
n simple zeros in K, ^ he a root-vector of f, 1 < p < oo and 0 < h < 1. 
Suppose x^ G IK” is an initial guess such that 

||a:° - ^llp < p = R{n,p,h)sep{f), (6.16) 

where R{n, p, h) is defined in fl6.13p . Then the Weierstrass iteration fll.3II is 
well-defined and converges guadratically to f with error estimates 06.141) and 

(6.17) 


Proof. It follows immediately from Theorem 16.51 and the obvious inequality 


E{x) < 


sep(/) 


which holds for every x G IK”. The estimate (I6.17p follows from the second 
estimate in (I6.14p and (I6.16p . □ 

Corollar vJ| 6 . 6 l is a generalization and improvement of the result of Dochev 
(see also j^). In the case p = cxo, he proved the estimate 06.17P under the 
initial condition 06.16p . 
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Corollary 6.7. Let f G ]K[z] be a polynomial of degree n > 2 which has 
n simple zeros in K, ^ 6e a root-vector of f, 1 < p < oo, 0 < h < 1 and 
0 < c < R{n,p) sep(/), where R{n,p) is defined in fl6.10p . Suppose G K” 
is an initial guess satisfying 


\\x°-^\\^<ch. (6.18) 

Then the Weierstrass iteration fll.Sp is well-defined and converges quadrati- 
cally to with error estimates fl6.14p and 

— flip < ch‘^’‘. (6.19) 

Proof. Obviously, the initial guess satisfies the initial condition (16.101) of 
Theorem 16.4[ Therefore, we obtain the estimates (16.lip which imply the 
estimates (16.141) . Indeed, we have 


A = (i){E{x^)) < (fi^ch) < h(f){c) < h(l){R{n,p)) = h 


since (p is quasi-homogeneous of the first degree on [0,1/2^^'^). The estimate 
(I6.19P follows from the second estimate in (I6.14p and (16.181) . □ 


Corollary |6.7| improves the result of Kyurkchiev and Markov 0 . In the 
case p = oo, they have proved the estimate (I6.19p under the initial condition 
(I6.16P but with a stronger condition for c. 


In 2002, Yakoubsohn 4l| published a y-theorem for Weierstrass method. 
He introduce the quantity 


7 (/) = max 7 (/,^i) where 7 (/,x) = max 

l<i<n k>l 


/(^)(x) 


k\f'{x) 


l/(fc-l) 


( 6 . 20 ) 


Recall that 'y{f,x) has been introduced by Smale in his famous work 35 


Corollary 6.8. Let f G ]K[z] be a polynomial of degree n > 2 which has n 
simple zeros in K, be a root-vector of f and 1 < p < oo. If & K” is an 
initial guess satisfying 


diO 


^ R{n,p, h) 

p~ 27 (/) 


( 6 . 21 ) 


where R{n,p,h) is defined in (I6.13p . then the Weierstrass iteration (II.3p is 
well-defined and converges quadratically to with error estimates (I6.14p . 
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Proof. It follows from Theorem 16.51 and the inequality 7(/) > 1/(2 sep(/)) 
which is due to Yakoubsohn [i^. □ 


Corollary 16.81 generalizes and improves the result of Yakoubsohu 4l| . In 


the case p = oo and h = 1/2, he proved the second estimate estimate in 06.141) 
under a stronger initial condition than 06.211) . 


7. Local convergence of the second kind of the Weierstrass method 


Let / G ]K[z] be a polynomial of degree n > 2 which splits in K, and let 
^ G K” be a root-vector of /. In this section we study the convergence of the 
Weierstrass method 01.31) with respect to the function of initial conditions 
E: V dehned by 


E{x) = 


X 


d{x) 


(7.1) 


where 1 < p < oo. Recall that P denotes the set of all vectors in K" with 
distinct components. We prove that the Weierstrass iteration function T 
dehned by 01 .Sp is also an iterated contraction with respect to this function 
E at the point The main result of this section generalizes, improves and 
complements the results of Wang and Zhao [^, Tilli and Han j^. 


Lemma 7.1. Let f G ]K[z] be a polynomial of degree n>2 which splits in 
K, ^ G K” be a root-vector of f, x & K” and 1 < p < oo. Then for i ^ j, 


\xi — > (1 — E{x)) di{x) and \xi — xf > dj{x), 

where E: V ^ M+ is defined by 07.ip . 

Proof. Setting in Proposition 15. If hi u = ^ and v = x and taking into account 
the dehnition of d(x), we obtain the hrst conclusion of the lemma. The second 
conclusion is obvious. □ 


Lemma 7.2. Let f G K[z] be a polynomial of degree n > 2 which splits in 
K, ^ G K” be a root-vector of f and 1 < p < oo. Suppose x G K” is a vector 
with distinct components such that 

fj{E{x)) > 0, (7.2) 

where the function E: T> M+ is defined by 07. ip and the real function is 
defined by 

= P'3) 
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Then f has only simple zeros in IK, Tx has pairwise distinct components, 


E{Tx)<p{E{x)) and \\Tx - ^ ^ f^iE{x))\\x - a, (7-4) 


where T\ T> C. K"" —)■ K"" is the Weierstrass iteration function defined by 
(na. and the real functions tp and fi are defined by ip{t) = t(3{t)/'fi{t) and 


m = 



t 

(n - ly/p 


n—1 

- 1 


(7.5) 


Proof. Note that condition fl7.2p implies that E{x) < since the func¬ 

tion is decreasing on M_|_ and ■^(1/2^/'^) < 0. Now it follow from Lemma 15731 
that the vector f has distinct components, which means that / has only sim¬ 
ple zeros in K. The second inequality in fl7.4p follows from fl6.7p taking into 
account that in this case \\u\\p < E{x). It remains to prove Tx eV and the 
hrst inequality in fl7.4p . Applying Proposition 15.21 with u = Tx and v = x, 
we get 


d{Tx) h {1- 2^/^ 


Tx — X 


d{x) 


d{x). 


(7,6) 


From the triangle inequality and the second inequality in fl7.4p . we obtain 


\\Tx - xii ^ \\Tx - eii + iix - eii A (1 + fi{E{x))\\x - eii 


which yields 

Tx — X 
d{x) 

Taking the p-norm, we get 

Tx — X 
d{x) 

Now it follows from fl7.6p that 


< E{x){l + fi{E{x)). 


d{Tx) ^ fi{E{x)) d{x). 


(7.7) 


It follows from this and fl7.2p that Tx eV. Combining fl7.7p with the second 
inequality in fl7.4p and taking into account fl7.2p . we get 


Tx-i 

d{Tx) 


P 0(i7(x)) 


d{x) 


(7.8) 


where the real function is dehned by (j){t) = I3{t)/fiit). Taking the p-norm 
in (17.81) . we obtain the hrst inequality in (I7.4p . □ 
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The next theorem is the hrst main result in this section. 


Theorem 7.3. Let f G ]K[z] be a polynomial of degree n > 2 which splits 
m K, ^ G K"" he a root-vector of f and 1 < p < oo. Suppose G K” is an 
initial guess with distinct components satisfying 

h{E{x^)) < 2, (7.9) 

where the function E is defined by fl7.ip and the real function h is defined by 

+ \ (7.10) 

Then f has only simple zeros in IK and the Weierstrass iteration lira is 
well-defined and converges to f with error estimates 

and ||x^-e|| (7.11) 

for all k > 0, where A = 4>{E{z^)), 6 = 'ijj{E{z^)), the real function is 
defined by (fit) = fi(t)/fi(t), and and {3 are defined by 07.31) and 07.51) . 
respectively. 

Moreover, if the inequality in 07.9p is strict, then the Weierstrass iteration 
converges quadratically to 

Proof. The function h is increasing and continuous on R+ with h(0) = 1 and 
/i(l/2^/^) > 2. Therefore, there exists a unique solution R = R{n,p) of the 
equation h{t) = 2 in the interval (0,1/2^/^h Hence, the initial condition 07.91) 
is equivalent to E{x^) G J, where J = [0, R]. It follows from h{R) = 2 that 

1 _ RID 

fi{R) = m) = Y:yfiR- (7-12) 

The function is decreasing on J with '0(0) = 1 and 4>{R) > 0. Hence, 

0 < '0(t) < 1 far all f G J (7-13) 

The function fi is increasing on J with 0(0) = 1 and 0(i?) < 1. From 
Example 12.21 we conclude that t0(t) is a strict gauge function of the second 
order on J. It is easy to see that (p is a quasi-homogeneous gauge function of 
the second degree on J. On the other hand, fl7.12p implies that i? is a hxed 
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point of (p. Then according to Proposition 12.41 (p is a gauge function of order 
r = 2 on J. 

It follows from fl7.13p that for every x with E[x) G J condition fl7.2p 
holds. Hence, Lemma [7.21 and Proposition 12.71 show that: 

• E\ T) ^ M_i_ is a functions of initial conditions of T with gauge function 
(p of order r = 2 on J. 

• T: P —)■ K” is an iterated contraction with respect to E at the point ^ 
with control function /3. 

• Every vector x^ E V satisfying condition fl7.9p is an initial point of T. 

Now the statement of the theorem follows from Theorem 13.31 □ 


In the next lemma we give a lower bound for the quantity R = R{n,p) 
dehned in the proof of Theorem 16.41 


Lemma 7.4. Let n > 2 and 1 < p < oo. Denote by R = R{n,p) the unique 
positive solution of the equation h{t) = 2, where the real function h is defined 
by (17.1 op . Then 


R > 


(n 


ly/q + ‘l^tq ■ 

This inequality becomes an equality if and only if n = 2 and p = 1. 


Proof. Let a = {n — b = 2fl^ and c = {n — 1)^/^. Note that a, b and c 
are greater than 1 and ac = n — 1. By the dehnition of R, we get 


(l + bR){l+R/c)"~' = 2. 


(7.15) 


Using Bernoulli’s inequality (1 + f)” > 1 + nt, we obtain 

f 1 + {a + b)R/n\^ / {a + b)R/n — R/c\^ ^ 1 + bR 

\ l + i7/c ) \ l + i?/c / “ l + i?/c 

From this and fl7.15p . we get 

{l + {a + b)R/n)'^ >2 (7.16) 

which proves (I7.14p . The equality in (I7.16p holds if and only if 6 c = 1 which 
is equivalent to n = 2 and p = 1. □ 

Theorem 16.41 together with Lemma YTM immediately implies the following 
result, which improves and complements the result of Han as well some 
other results. 


25 
































Corollary 7.5. Let f G ]K[z] be a polynomial of degree n > 2 which splits 
m K, ^ G K" he a root-vector of f and 1 < p < oo. If E K” is an initial 
guess with distinct components satisfying 


1)1/, + 2V. ■ 

then f has only simple zeros in K and the Weierstrass iteration mi is 
well-defined and converges to f with error estimates (I7.1ip . Moreover, the 
convergence is quadratic provided that n > 3 or p > 1 or fl7.17p holds with 
strict inequality. 

Corollary 7.6. Let f G K[z] be a polynomial of degree n >2 which splits in 
K, G K” be a root-vector of f and 1 < p < oo. If E K"’ is a vector with 
distinct components satisfying 


E{x°) = 


X 


d{x^ 


E(x°) 


d{x^) 


< R 

p 


2l/(n-l) _ ^ 


(7.18) 


then f has only simple zeros in K and the Weierstrass iteration fll.Sp is 
well-defined and converges quadratically to f with error estimates (17.111) . 

Proof. It follows from Corollary 17.51 because the right side of fl7.17p is less 
than the right side of fl7.18l) □ 

Corollaiy (17.6p is a generalization and improvement of the result of Wang 
and Zhao [37|. They have proved the convergence of the Weierstrass iteration 
fll.Sp for a polynomial f E C[z] under the stronger initial condition 


X 


V2- 1 

4 " V2 - 3 


min di(x^). 

l<i<n ^ ’ 


than fl7.18p in the case p = oo . 

The next theorem is the second main result in this section. 


(7,19) 


Theorem 7.7. Let f G K[5;] be a polynomial of degree n>2 which splits in 
K, G K” be a root-vector of f and 1 < p < oo. Suppose c: [0,1) —)■ M+ is 
a nondecreasing function such that 

c{t) <{n-l)RP{{l + t)R’^^-^^-1) for all tG[0,l). (7.20) 
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Suppose also there exists a G (0,1) sueh that 


t c(t) < c(t^) (1 — 2^/'^ (1 + t) c(t) for all tG[0,(T], (7.21) 


where 1 < (? < cxd and 1/p + 1/q = 1. Let G K” be a vector with 
components satisfying 


X 


d{x^ 


< c(a) 


Then f has only simple zeros in IK and the Weierstrass iteration 
well-defined and converges quadratically to f with error estimates 


distinct 

(7.22) 
(|1.3p is 


-^W^a^'Wx'^ -^\\ and ^ (7.23) 


for all k > 0. 

Proof. For the sake of simplicity, let b = 2^/^. It is easy to show that 
1 — 6 (1 + t) c{t) > 0 for all t G [0,1). Define the functions fj and fd by 07.31) 
and 07.5p on the interval J = [0,c(cr)]. It is easy to show that the function 
^|J can also be rewritten in the form 

fjit) = 1- 2^^H{1 + ^{t)). (7.24) 

Condition 07.201) can be rewritten in the following equivalent form 

/9(c(t)) < t for all t G [0,1). (7.25) 

Note that (d is an increasing function on J with values in [0,1) since 

jdit) < /3(c((t)) < a < 1 for all t G J. 

The function '0 is a decreasing function on J with values in (0,1] since 

> '0(c(c’‘)) = l—bc{a){l-\-(d{c{a))) > 1—b (1+cr) c{a) > 0 for all f G J. 

Therefore, we can define the function 0 on J by (p{t) = t fd{t)/'ijj{t). Condition 
07.2ip implies the following one 

7’(c(t)) < c(t^) for all t G [0,cr]. (7.26) 

Indeed, it follows from 07.24p . 07.251) and 07.211) that 


p{c{t) 


c{t)P{c{t) tcjt) 2 

i/{c{t)) ~ l-b{l-\-t)c{t) ~ 
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Note also that the function if is nondecreasing on J and f{J) C J since 

< <d(c(o')) = c(cr^) < c{a) for all f G J. 

Consider again the Weierstrass iteration function T: D —)■ K"' and the func¬ 
tion E: V ^ M+ dehned by fIT.ip . For every x E D with E{x) G J condition 
07.21) holds. Indeed, we have 

Tjj{E{x)) > ^p{c{a)) > 0 . 

It follows from Lemma 1712] that E is a. functions of initial conditions of T with 
gauge function f on J, and that T is an iterated contraction with respect to 
E at the point ^ with control function f3. Now the statements of the theorem 
follow from Theorem 13.51 □ 


There are a lot of functions that satisfy condition fl7.20p of Theorem 17.71 
For example, each of the functions 


c{t) = 


2t — r 


and 


c{t) = 


2t 


2 (n-1)1/9 (n - 1)1/'? (f + 2) ■ 

satishes (I7.20p . This statement follows from the obvious inequalities 

2t 


(7.27) 


(1 t) 




1 > and ln(l + t) > 


n 


t+2-^ Y' 


Applying Theorem 17.71 with the hrst function c: [0,1) M+ dehned by 
fl7.27p we obtain at the following result. 

Corollary 7.8. Let f G IK)^;] be a polynomial of degree n > 2'?+i -|- 1 which 
splits oner K, ^ be a root-vector of f, 1 <p < oo, and let a G (0,1/2] he such 
that 


CT 


l)(2-a)(2 


a 


1 — cr 


< 2 


n 


1/9 


Suppose x^ G K" is a vector with distinct components satisfying 


x^ -f 


d{x^ 


< 


2a — 

2 (n- 1)1/9 


(7.28) 


(7.29) 


Then f has only simple zeros in IK and the Weierstrass iteration (II.3p is 
well-defined and converges guadratically to the vector ^ with error estimates 


X 


k+l 


eii ^ 




^11 and 




ell, (7.30) 
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11^^ ~ ^lloo ^ 0’^*' max di{x^) (7-31) 

l<i<n 

for all k>0. Besides, if n > + then condition fl7.28l) can be dropped. 

Proof. The statement follows from Theorem 17.71 with the second function 
c: [0,1) —)■ M+ dehned by fl7.27p . The estimate fl7.3ip follows from the sec¬ 
ond estimate in fl7.30p and the initial condition fl7.29p . It should be noted 
only that condition fl7.28p is equivalent to fl7.2ip and that fl7.28p is satished 
automatically if u > 2^^+^ -|- 1. □ 


Corollary fl7.8p is a generalization and improvement of the result of Tilli 
3^ who has proved the error estimate fl7.3ip under the condition fl7.29p with 


p = cxo for a polynomial f E C[z] of degree n > 9. 


Applying Theorem 17.71 with the second function c: [0,1) —?■ R+ dehned 
by fl7.27p we arrive at the following result. 


Corollary 7.9. Let f G IK)^;] be a polynomial of degree n > 2*?+^ -|- 1 which 
splits over f he a root-vector of f and 1 < p < oo. Suppose a is a real 
number satisfying 

^ ^ ^ in- - 21 +iA 

^ ^ ^ - {n- 1)1/9 + 21+V9 

and G K” is a vector with distinct components such that 


d(x°) 


2a 

(n — 1)^/1 (a -I- 2) 


(7.32) 


Then f has only simple zeros in IK and the Weierstrass iteration fll.Op is 
well-defined and converges guadratically to f with error estimates fl7.30l) . 


Corollary 17.81 can be state in the following equivalent form. 


Corollary 7.10. Let f G IKfz] be a polynomial of degree n > 2^^^ + 1 which 
splits over K, ^ he a root-vector of f and 1 < p < oo. Suppose G K"" is a 
vector with distinct components satisfying 


X 


d{x° 


< 


2 [{n-iy/i - 2 ^+ 1 /'?] 

(n - 1)1/9 [3(n - 1)1/9 + 21+1/9] ‘ 


Then f has only simple zeros in IK and the Weierstrass iteration 
well-defined and converges guadratically to f with error estimates 


\x‘ 


k+l 


-ell ^ A" ||x^-e|| and llx'^-eil ^ A" “ix^'-eil 


(7.33) 
01.3p is 

(7.34) 
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for all k > 0, where 


A = 


2(n- 1)1/3 E(a;°) 

2 - {n - E{xO) 


8. Semilocal convergence of the Weierstrass method 

In this section, we prove a new convergence theorem for the Weierstrass 
method under computationally verihable initial conditions. The main result 
of this section generalizes, improves and complements all previous results 

and Zhao 
Petkovic, 


381 |44| . Petkovic, Carstensen and Trajkovic [l_^, 

Herceg and Ilic 211, Batra j^, Han [^, Petkovic and Herceg j20| and Proinov 
24| . The new result in this section also gives computationally verihable error 


in this area, which are due to Presic 


Zheng l45l 


Wanj 
Petkovic llj 


estimates, a localization formula for the Weierstrass iteration (OD as well 
as a sufficient condition for a polynomial / G K[z] of degree n > 2 to have n 
simple zeros in the held K. Finally, we provide an example which shows the 
exactness of the error estimates of our semilocal theorem for the Weierstrass 
iterative method. 

We study the convergence of the Weierstrass method (II.dh for a polyno¬ 
mial / G ]K[z] with respect to the function of initial conditions E: V ^ M+ 
dehned by 

W{x) 


E{x) = 


d{x) 


{1 < p < oo) 


( 8 , 1 ) 


We prove that the Weierstrass iterative function T: T) C. K"" —)■ K” dehned 
by fll.Sp is an iterated contraction with respect to E. 

We begin this section with a well-known result whose proof we include 
for completeness. 

Proposition 8.1. Let f G IK)^;] be a monic polynomial of degree n >2, and 
let X G K” be a vector with distinct components. Then for all 2 ; G K, 

n n 

/(^)=n 

i=l j^i j=l 

Proof. Applying Lagrange’s interpolation formula to the polynomial 

n 

= /(^) - - Xj) 

i=i 

at the nodes Xi,..., x^ we get the desired presentation of /. □ 
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Using the Weierstrass correction Wf one can state the following basic 
existence result for polynomial zeros. 

Proposition 8.2 (Basic Existence Theorem). Let f G be a polynomial 
of degree n > 2, and let {x^) be an infinite sequence of vectors in K” with 
distinct components satisfying the following two conditions: 

(i) The sequence (x^) converges to a point f in K”; 

(ii) The sequence {Wf{x^)) converges to the zero-vector in K"'. 

Then f is a root-vector of f. 

Proof. Without lose of generality we may assume that / is a monic polyno¬ 
mial. Applying Proposition 18.11 with x = we obtain 

n n 

m =n +n ("= - 

i=l j^i j=l 


for all k. Passing to the limit when /c —)■ cxo, we get f{z) = YYj=i — ij) 
which completes the proof. □ 

The following remarkable result for the Weierstrass iteration is an imme¬ 
diate consequence of Proposition 18.21 


Proposition 8.3. Let f G be a polynomial of degree n > 2. If for 

some initial guess x° G K”, the Weierstrass iteration fll.3p is well-defined 
and converges to a vector f G K”, then f is a root-vector of f. 


Let us g ive some historical historical notes about Proposition [8]3l In 1972, 


Petkov |l6l. p. 272] has briefly noted (without proof) that if the Weierstrass 
iteration (II.3p for a complex polynomial / converges to a vector ^ G C”, then 
^ is a root-vector of /. This result for complex polynomials with simple zeros 
has been proved by Zheng in 1982 and re-obtained in 1994 by Hopkins, 
Marshall, Schmidt and Zlobec jsi Proposition 4.1]. 


Proposition 8.4. Let $: 77 C K"' —)■ K” be a mapping defined on a set D 
which contains only vectors with distinct components, and let E: 77 —)■ M_|_ be 
defined by 


E{x) 


<h(x) 

d{x) 


(1 < p < c)o). 
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Let X ^ D be such that E{x) G J, where J G ^ is an interval containing zero. 
Furthermore, let there exist two function J ^ M+ such that 

I3{t) < 1 — 2^^‘^t'y{t) for all t e J. (8.2) 

Then the closed disks 


Di = {z E K : \z — Xi\ < r,}, i = 1, 2,..., n, 


(8.3) 


where 


Ti = 


l{E{x)) 

l-P{E{x)) 




are mutually disjoint. 

Proof. For simplicity, we set b = 2}!^ and C = 7 (i?(x))/(l — ft{E{x))). It 
follows from fl8.2p that 


bt-f{t) 

1 - m 


for all t E J. 


This implies that bC E{x) < 1 since E{x) E J. To prove that the disks (18.31) 
are mutually disjoint it is sufficient to show that Vi + Vj < \xi — xf, that is, 


C(|$j(a;)| + |$j(a;)|) < -Xjl for i^j. (8.4) 


Suppose C > 0 since the case C = 0 is obvious. From the dehnition of d{x), 
the power mean inequality Mi < Mp and the inequality bC E{x) < 1, we 
obtain 


$^(x)| + |<l>j(x)| ^ |<F^(a))| ^ |4>j(x)| 
\xi — Xj\ ~ di{x) dj{x) 


< bE{x) < 


1 

C 


which yields (I8IT|) . 


□ 


Lemma 8.5. Let f E IK)^;] be a polynomial of degree n>2 which has n simple 
zeros inK, ^ E IK” be a root-vector of f, x E K” and 1 < p < oo. Then 


Xi — > (1 — 2^^'^E{x)) dj{^) and \xi — > (1 — E{x)) di{^) 


for i ^ j, where F': K” —)■ M+ is defined by (16.11) . 

Proof. Setting in Proposition 15.11 u = x and v = ^ and taking into account 
the dehnition of d{^), we obtain the statement of the lemma. □ 
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Lemma 8.6. Let f G ]K[z] be a polynomial of degree n>2 and 1 < p < oo. 
Suppose X G K"" is a vector with distinct components satisfying 


W{x) ^ 


d{x) ^ * 


(8.5) 


Then the vector x = x — W{x) has distinct components and 



= ( 8 . 6 ) 


Proof. Applying Proposition 15.11 with u = x and v = x and taking into 
account condition fl8.5p . we conclude that x* ^ Xj and Xj ^ Xj for i ^ j. 
Hence, both sides of fl 8 . 6 p are well-dehned. In particular, the vector x has 
distinct components. Let i G be hxed. If Xj = Xj, then fl 8 . 6 p holds trivially. 
Assume that Xj 7 ^ Xj. It follows from Proposition 18.11 that for every G IK 
such that z ^ Xj (j = 1 ,..., n). 



which can be rewritten in the form 



Setting here z = Xi and taking into account that Xi = Xi — IPj(x), we get 



Therefore, 



which coincides with fl 8 . 6 p . 


□ 
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Lemma 8.7. Let f G ]K[z] be a polynomial of degree n>2 and 1 < p < oo. 
Suppose X G K"" is a vector with distinct components such that 


E{x) = 


W{x) 


d{x) 


< 


2V9 


Then the vector Tx = x — W{x) has distinct components and 

E{Tx) < g}{E{x)) and \\Tx — T‘^x\\ ^ (3{E{x)) \\x — Tx\\, 
where the real functions (p and ft are defined by 


p{t) = 


{n — 

(1 -t)(l - 2Vn) 
(n — t 


1 + 


(n — ly/p (1 — 2^/*? t) 




i-t 


1 + 


(n - ly/p (1 - 2Vn) 


n—1 


n—1 


Proof. Proposition IS.lT L with u = Tx and v = x yields 

d{Tx) ^ if{E{x)) d{x), 
where the real function ip is dehned by 


(8.7) 

(8.8) 

(8.9) 

( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 


Obviously, (18.1 1 1I implies that Tx has distinct components. For the sake of 

simplicity, we use the following notations: 

x = Tx, a={n — iy^‘^, 6 = 2^/'^ and c=(n — 

It follows from Lemma 18.61 that 

||lF(ri)|| ^ffp||lF(K)|l, (8.13) 

where 

Wj(x) 


a = max 
iGL. 


E 


OC 2 j 


and Li = max 11(1 

iG/n J-J. V 




wyx) 


X 7 X 7 


Applying Proposition 15.11 with u = x and v = x, we get 

\xi — Xj\ > (1 — E{x)) dj{x) and \xi — Xj\ > (1 — bE{x)) dj{x). 
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Using the power mean inequality Mi < Mp and Proposition 15.51 we obtain 


a < 


aE{x) 


and /i < ( 1 + 


E{x) 


c (1 — bE{x)) 


1 — E{x) 

Now from 08.131) and 08.14p . we get 

\\W{Tx)\\^P{E{x))\\W{x)\\ 


n—1 


(8.14) 


(8.15) 


which coincides with the second inequality in 08. 8 p since W{x) = x — Tx. It 
follows from 08.151) and 08. lip that 


W{Tx) 

. m{x)) 

W{x) 

d{Tx) 

- y{E{x)) 

d{x) 


Taking the p-norm here, we get the hrst inequality in 08.8p since (3 = ^ 

Now we can state and prove our semilocal convergence theorem for the 
Weierstrass method which is the main result of this section. 


Theorem 8.8. Let K be a complete normed field, f G ]K[ 2 ;] be a polynomial 
of degree n > 2 and 1 < p < oo. Suppose x^ G K" is an initial guess with 
distinct components satisfying 

E{x^) < 1/2^/" and (p{E{x^)) < 1, (8.16) 


where the function E is defined by 08.ip and the real function is defined by 


(p{t) = 


(n 


Ifi/H 


1 + 


(1 - t)(l - 2^11) V (n - ly/P (1 - 2^/n) 


n—1 


(8.17) 


Then the following statements hold true. 

(i) Convergence. Starting from x^, the Weierstrass iteration 01.3p is 
well-defined, remains in the closed ball U{xo,p) and converges to a root- 
vector f of f, where 

\\W{xy\\ 

^ l-fi{E{xy) 

and the real function fi is defined by 08.10p . Besides, the convergence 
is quadratic provided that 0(U'(a:°)) < 1. 
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11 


A PRIORI ESTIMATE. For all n > 0 we have the estimate 


nk \2^-l 


(8.18) 


111 


where A = (j){E{x^)), 9 = 'ijj{E{x^)) and the real functions fj is defined 
by dHUD. 

First a posteriori estimate. For all k > 0 we have the following 
estimate 

lk"-eil A 


I 


IV 


l-P{E{x^))' 

Second a posteriori estimate. For all k >0 we have the following 

^fcAfc 


estimate 


\X' 


kfi-1 


eii A 


I I 


l-OkiXkfi 

where Xk = (j){E{x^)), 9k = 'if{E{x’^)). 

(v) Some other estimates. For all k >0 we have 


( 8 . 20 ) 


xk+2 _ ^fc+l|| ^ 


k+1 


VI 


( 8 . 21 ) 
( 8 . 22 ) 

Localization of the zeros. If (j){E{x^)) < 1, then f has n simple 
zeros in IK. Moreover, for every k > 0 the closed disks 


X 




\x — X 


Df = {zeK:\z-x’:\<rf}, * = 1, 2,..., n, (8.23) 


where 


,.k ^ |M^i(T)| 


1 - ;3(B(i‘) ’ 

are mutually disjoint and each of them contains exactly one zero of f. 


Proof. Denote hj R = R{n,p) the unique solution of the equation (j){t) = 1 in 
the interval (0, 1 / 2 ^A^_ Then the initial conditions (I8.16p can be rewritten in 
the form E{x^) < R. The function ip dehned by fl8.9p is quasi-homogeneous 
of the second degree on [0,1/2^/'?) and R is its hxed point. It follows from 
Proposition 12.41 that <yC is a gauge function of the second order on J = [0, A]. 
The function fi is increasing on J satisfying (3 = and (3{R) = fi{R) < 1- 
Hence, the function t(3(t) is a strict gauge function of the second order on 
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J. Now we shall apply Corollary 14.41 to the Weierstrass iteration function 
T \ V C. K"" —)■ K"'. From Lemma [8.71 and Proposition 12.71 we conclude that: 

• E\ T) ^ M+ is a functions of initial conditions of T with gauge function 
ip of order r = 2 on J. 

• T: P —)■ K"’ is an iterated contraction with respect to E with control 
function j3. 

• Every vector satisfying condition fl8.16p is an initial point of T. 

Now it follows Corollary 14.41 that the Weierstrass iteration fll.Sp is well- 

dehned and converges to a vector ^ G K”. The claims (i)-(v) follow im¬ 
mediately from Corollary 14.41 and Proposition 18.31 It remains to prove the 
claim (vi). Let E{x^) < 1/2 and 0(i7(x°)) < 1, that is, E{x^) < R. Then 
E{x^) < R since every x^ is an initial point of T. Taking into account that 
/3(f) = and < 1 on [0, R), we conclude that /3(f) < '0(f) for all 

f G [0,i?). Therefore condition fl8.2p holds with 7 (f) = 1. Applying Proposi¬ 
tion |8]1] with <F = IT we conclude that the disks fl8.23p are mutually disjoint. 
On the other hand, it follows from claim (iii) that each of these disks contains 
at least one zero of /. Therefore, each of the disks fl8.23p contains exactly 
one zero of /. In particular, / has n simple zeros in K. This complete the 
proof of the theorem. □ 


It should be noted that recently Proinov and Petkova |33(] proved another 
semilocal convergence theorem for the Weierstrass method under initial con¬ 
ditions involving the Viete operator. 


Remark 8.9. Let R{n,p) be the unique solution of the equation 0(f) = 1 in 
the interval (0,1/2^/'^), where 0 is dehned by (I8.17p . Then the initial condi¬ 
tions fl8.16p of Theorem 18.81 can also be written in the following equivalent 
form 

E{x^) < R{n,p), (8.24) 

where E\ T> ^ M+ is dehned by (18.ip . According to Theorem l 8 . 8 l the domain 
of quadratic convergence of the Weierstrass method is the set 


= {x eV ■. E{x) < R{n,p)}. (8.25) 


Below we show that this cannot be improved. 
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Example 8.10. In this example we show that if p = oo, then the domain 
of quadratic convergence 08.251) of the Weierstrass method as well as the 
estimates OS.lSp . 08.19p . O8.20p . 08.2ip . fl8.22p cannot be improved in the 
sense that for every polynomial f{z) = + aiz + 02 in ]K[z] with zero 

discriminant, there exist inhnitely many initial guesses G such that 
E{x^) = R{2, 00 ) = I and the following two statements hold true: 

(a) The Weierstrass iteration 01.3p is well-dehned and converges linearly to 
the root-vector ^ = (a, a) of /, where a = —ai/(2ao). 

(b) Each of the estimates 08.18p - 08.22p becomes equality. 

To prove this, we choose a vector = {a + b,a — b), where b is an arbi¬ 
trary element of K. It is easy to calculate that E{x^) = R{2, 00 ) = where 
E is dehned by 08.ip with p = 00 . According to Theorem 18.81 the Weier¬ 
strass iteration 01.3p is well-dehned and converges to a root-vector of / 
with estimates 08.18p - 08.22p . which can be written in the following form: 


xk+2 _ ^fc+l|| 
imfc+i _ ^fcii 



(8.26) 


On the other hand, the Weierstrass iteration 01.3p for f{z) = qqz'^ + aiz + a 2 
with the initial guess = {a + b,a — b) can be written in the form 


X 


k 



Hence, the sequence [x^) converges to the vector ^ = {a, a). Besides, 

f = i (x^ - {). 

This shows that the order of convergence {x^) is exactly one (with asymptotic 
constant which proves the statement (a). Furthermore, it is easy to check 


38 




























that {x^) satisfies the following identities: 


„k+l 


x^ -e 


x^ 

= 2(x'=-x^+^), 

x^ -e 

= x'^-'-x'=). 

- x^+^ 

= ^(x'^-x'^+i). 


= (i)‘ 


These identities show that all the inequalities fl8.26p become equalities. This 
proves the statement (b). □ 


Remark 8.11. It follows from Theorem 18.81 and Lemma ISTl that if an initial 
guess G K” satisfies the initial conditions fl8.16p for some 1 < p < oo, 
then the Weierstrass iteration fll.3p is well-defined and can be presented in 
the following two-point form: 

x^ = x^ — Wf{x^), 

_ rfk rr.k _ ryk-l 

=4 - (x‘ - xp‘) 5: n (*=1..... n) 

^3 j^i -^3 

(8.27) 

k = 1 , 2 ,... 

Let us note that in the Weierstrass method fll.Sp . we need to compute the 
polynomial values /(xf) for each iterative step. Under the initial conditions 
fl8.16p the two-point iterative method fl8.27p is equivalent to the Weierstrass 
method fll.3p but no longer needs to evaluate /(xf) for each iteration after 
the first step. The two-step iterative process 08.271) was first presented in 
1964 by Dochev and Byrnev in a slightly different form (see also Zheng 
451 and Yao 42 ). 


Remark 8.12. In an earlier work 2^ (see also 25|]) we have stated without 


proof a weaker version of Theorem 18.81 In 2^ we give a detailed comparison 
of our old theorem with previous results. It should be noted that all corol¬ 
laries given in 2^ can be improved using Theorem 18.81 instead of Theorem 1 
of 2^. We end this section with a result which improves Corollary 2 o f 124 . 
This result generalizes and improves the results of 0 M, 0, m- 
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Corollary 8.13. LetK be a complete normed field, f G ]K[z] be a polynomial 
of degree n > 2 and 1 < p < oo. Suppose G K"" is an initial guess with 
distinct components satisfying 


W{x^) 

d{x^) 


2{n 


1 

1)V? + 2‘ 


Then all conclusions (i)-(v) of Theorem {KM hold true. Moreover, ifn>3, 
then f has n simple zeros in K, the Weierstrass iteration ma converges 
quadratically to f and for every k >0 the closed disks fl8.23p are mutually 
disjoint and each of them contains exactly one zero of f. 
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